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ABSTRACT 

Superbursts are thought to be powered by the unstable ignition of a carbon-enriched layer formed from the 
burning of accreted hydrogen and helium. As shown by Cumming & Bildsten, the short recurrence time hinges 
on the crust being sufficiently hot at densities p > 10 9 gcm~ 3 . In this Letter, we self-consistently solve for the 
flux coming from the deep crust and core. The temperature where the carbon unstably ignites is only weakly 
sensitive to the composition of the ashes of H/He burning, but does depend on the thermal conductivity of the 
inner crust and the neutrino emissivity of the core. The observed superburst recurrence times and energetics 
suggest that the crust thermal conductivity is low, as if the crust were amorphous instead of crystalline. If the 
conductivity is higher, such as from a lattice with impurities, then matching the superburst properties require 
that the neutrino emissivity be not stronger than modified Urea. Observations of superbursts — energetics, 
recurrence times, and cooling times — therefore complement observations of isolated cooling neutron stars and 
soft X-ray transients in constraining properties of dense matter. Perhaps the most interesting object in this 
regard is KS 173 1 —260, which produced a superburst during its protracted accretion outburst but had a rapidly 
declining quiescent luminosity. 

Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars: neutron — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

With the launch of BeppoSAX and RXTE, long-term moni- 
toring of low-mass X-ray binaries became feasible. One fruit 
of this was the discovery of superbursts: these have a rapid 
rise and long decay, analogous to type I X-ray bursts, but 



are roughly a thousand times more energetic i Corndisse_etal 



times more energetic ( 
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2000; in'tZandetal 
2002; Kuulkersl l2002 : 
2001). The recurrence time is of order a year; for one 
source, 4U 1636 —53, two superb ursts were observed sep- 
arated by 4.7 yr (Wiinands 2001). Although most super- 
bursts are observed at mass accretion rates 0.1-0.3 M^dd, 
Mgdd ~ lO^gs -1 being the Eddington accret ion rate, re- 
cently in't Zand, Cornelisse, & Cumming (2004) detected su- 
perbursts from the rapidly accreting source GX 17+2. 

Despite their similarities to type I X-ray bursts, the 
fuel for the superburs ts was not immediately apparent. 
Taam & Picklumi ( f!97l had originally investigated the ther- 
mal stability of a pure 12 C layer and found that such a 
layer would unstably ignite when the accumulated mass 
was ~ 10 27 g. This produces a very energetic burst, E ~ 
10 44 ergs, with a r ecurre nce time of 300 yr (M/O.lMedd)- 
lBrown & Bildstenl |l998) demonstrated that carbon flashes 
with greatly reduced recurrence times could occur on a neu- 
tron star (NS) accreting pure He at a locally super-Eddington 
rate, such as might happen on an accreting X-ray pulsar. This 
explanation did not fit the observed superbursts because none 
of the sources showed pulsations in the persistent emission 
and, with the exception of 4U 1820—30, all were thought to 
accrete a H/He mixture t hat burned to matter with A ~ 100 
fcchatz et all 1 9991 l200l. 

Cumming & Bildsten (2001) suggested a resolution to this 
puzzle. They noted that only a small 12 C mass fraction 
X( 12 C) ~ 0.1 is sufficient to produce a thermonuclear run- 
away, and that the reduced conductivity of a heavy element 
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mixture (e.g., 0.9:0.1 104 Ru: 12 C) would produce a steeper 
thermal gradient for a fixed flux emergent from the deeper 
crust. The flux from the deep crust is powered by elec- 
tron captures, neutron emission s and pycnonuclear reactions 
( lHaensel & Zdunikll990tl2003l) that release roughly 1 .4 MeV 
per accreted nucleon; the outward-directed luminosity from 
the deep core is roughly L nuc < 0.1 MeV (M/m„), for steadily 
accret ing NSs with modified Urea neutrino cooling in their 
cores (Brown 2000). Here m u is the mean nucleon mass. Un- 
der these conditions, superbursts would ignite for an accumu- 
lated mass ~ 10 25 g, which is consistent with the observed en- 
ergetics and recurrence times. The superbursts observed from 
GX 1 7+2 are consistent with this scenario din't Zand et alJ 
l200l . 

In this Letter, we calculate the thermal structure of the 
entire crust using a range of core neutrino emissivities and 
crust thermal conductivities. The flux in the carbon-enriched 
layer is not fixed, but increases with (A)/(Z 2 ), where (A) 
and (Z 2 ) are the mean nuclear mass and squared nuclear 
charge of the carbon-enriched layer. As a result, the tem- 
perature at which superbursts ignite is relatively insensitive 
to the dominant isotope formed during H/He burning. Fur- 
thermore, the temperature of this layer is sensitive to the na- 
ture of the thermal conductivity and neutrino emissivity of the 
deep crust and core. As a result, the superburst properties — 
recurrence times, energies, and cooling times — open a new 
view into the interior physics of NSs, complementary to ob- 
servation s of isolated cooling neutr on stars (for a recent re- 
view, see |\akovkiv_&Pethic^2004j ) and soft X-ray transients 
llRutledge et alJ 1200 It lYakovlev et aU 120031: [jakovlev etalJ 
2004). In particular, the obser ved energetics and recurrence 
times agree with the scenario of Cummi ng & Bildstenl ( 120011) 
only if the nuclei in the NS crust are disordered or amorphous. 

2. THE THERMAL STRUCTURE OF THE CRUST 

To calculate the c rust th ermal structure, we use the method 
outlined in Brown (2000) with a few modifications. First, 
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we take the core neutrino emissivity as a free parameter, 
rather than compute it self-consistently for a given equation 
of state (EOS). This is because our goal is to study the crust 
thermal structure for a given average core neutrino emissiv- 
ity . To simplify our ca lculations, we follow the approach 
of lYakovlev & HaenseD d2003l> and use an analytic al den- 
sity prescription p = p D [l — (r/R) 2 } dTormanlll939l) . where 

p = 15{87t)- l MR- 3 with R = (2GM/c 2 ) [l - (1 +z) _2 j 

and l+z = [1 — 2GM / (Rc 2 )]~ l l 2 . This prescription gives 
a reasonable approximati on to the p(r) foun d using modern 
nuclear EOS's (see Lattimer & Prakash 2001, and references 
therein). With this choice for the EOS, the NS's mass M 
and redshift z specify its mechanical structure. We locate the 
crust-core boundary by a Maxwell construction using an EOS 
for the inner crust {N egele & Vautheri nll 19731) and the core 
( Akmal, Pandharipan dT&Ravenhalll ll998). which gives the 
transition density Pb = 1.6 x 10 14 gcm~ 3 . 

With the core mass and radius specified, we compute the 
structure of the crust by integrating the general relativis- 
tic equation of hydrostatic balance from P(p = Pb) to P = 
g x 10 9 gcm~ 2 . We then solve the thermal equations 
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by a relaxation method JPress et alJ 119921) . Here T and 
L denote proper values and e v « (1 +z)~ 2 - The condi- 
tion at the crust/core boundary is L(p = pt) + L v =0: the 
heat conducted into the core must be re-radiated as neutri- 
nos. At a column depth y = 10 9 gcm~ 2 we set T — 2.5 x 
10 8 K, in accordance with es timates for mixed H/He ignition 
JCumming & Bildstenl2000l) . 

The crust thermal profile is controlled by the EOS, the 
volumetric neutrino emissivity e v and nuclear heating rate 
£ nU c, and the thermal conductivity K. The crust EOS 
has contributions from degenerate, relativistic electrons, for 
which we use a tabulati on of the Helmholtz free energy 1 
( Tim mes & S we stvl 120001): from ions, for which we use the 
free-energy fits of Farouki & Hamaguchi (1993); and from 
free ne utrons, for which we us e a compressible liquid-drop 
model (M ackie & Bavml| 19771) . The neutrino emission in 
the inner crust is dominated by neutrino pair bremsstrahlung 
(Maxwell 1979 j), for which we use t he rate calculated by 
lHaensel. Kaminker. & Yakovlevl in'996). As in Br ownl (120001) . 
we do not resolve the reaction layers in the deep crust, but 
instead use a distributed e nuc normalized so that the total lu- 
minosity from the deep crust reactions is L nuc = (1.4MeV) x 
(M/m u ). 

Heat is conducted by electrons, and the thermal conductiv- 
ity is given by the Wiedemann-Franz law, with the conduc- 
tivity proportional to an electron-ion relaxation tim e T, fo r 
which we use the fitting formula of Potekhin et al. ( 1999), 
in the mean-ion approximation. A substantial uncertainty in 
computing the transport coefficients of the crust is its compo- 
sition. During the lifetime of a low-mass binary, it is possi- 
ble to replace a substantial portion of the ~ 0.01 M Q crust. 
The accreted hydrogen an d helium burn via th e rp-process 
(Walla ce & Woosleviri981l) : both st eady-state (ISchatz et all 
1999) and time-dependent, one-zone (Schatz et al.l200ll) cal- 
culations predict that the ashes are a broad distribution 

1 Extension to nuclear densities courtesy F. Timmes 



of ions peaking with (A) w 100. Recent calculations of 
one-dimensional models of episodic accretion and unsta- 
ble H/He ignition for a range of accretion rates and metal- 
licities find that repeated bursts drive the me a n weight of 
the a shes to (A) sa 64 dWooslev et alJ 120041 iFisker et ail 
1 20041). For definiteness, w e use the compositional profile of 
Haensel & Zdunik (2003), which presumes that the ashes of 
H/He burning have (A) « 100; in t he inner crust the co mposi- 
tion is similar, however, to that of Haensel & Zdunik ( 1990), 
which assumed the a shes of H /He bur ning are 56 Fe . 

As pointed out by Cumming & Bildsten (2001), the com- 
position of the carbon-enriched layer can have a strong effect 
on the temperature profile. In this layer the ions are not crys- 
tallized and the conductivity scales roug hly as (A)/(Z 2 ). As 
a result, one would expect a steeper gradient for a heavier ion 
mixture and a given flux. The flux is not fixed, however, and 
must be determined self-consistently. Figure ^ shows the lu- 
minosity emergent from the crust for different values of Z 2 /A 
in the layer 10 9 gcm~ 3 < y < 10 13 gcm~ 3 . The connected 
squares mark the results of the calculations; the solid line in- 
dicates the slope A/Z 2 . The luminosity is scaled by a fidu- 
cial value, L c = 1 MeV x (M/m u ) = 9.6 x lO^ergss" 1 Mi 7 , 
and the numbers beside each line indicate the accretion rate 
Mn =M/(10 17 gs- 1 ). 

To understand why the thermal gradient depends weakly on 
the composition of the layer below the H/He burning shell, 
consider a simple two-zone model of the crust, with the heat- 
ing from the pycnonuclear reactions at the interface between 
two zones. When most of the heat generated in the crust 
flows into the core, the maximum temperature in the crust, 
and therefore the thermal gradient between the H/He burning 
layer and the electron capture layers, is roughly fixed. This 
approximate constancy in the temperature gradient implies 
that the combination L(Z 2 ) /(A) is roughly constant (Fig. 0. 
The superburst ignition temperature depends weakly on the 
composition of the carbon-enriched layer. At lower accre- 
tion rates and smaller Z 2 /A, this simple argument weakens, 
as the luminosity in the outer layers becomes a larger frac- 
tion of Lnuc- For use in future computations, we fit the lu- 
minosity in the carbon-enriched layer as a function of M\i 
and Z 2 /A for 1 < M 17 < 10 and 8 < Z 2 /A < 20. We write 
L = qMN A (lMeV) = (9.65 x lO^ergss" 1 ) x qM l7 , with 



q = 0.386 
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This fit has an rms error of 1.7%, with a maximum error of 
5.1% for .21% = 10 and Z 2 /A = 20. 

3. SENSITIVITY TO THE PHYSICS OF THE DEEP CRUST AND 
CORE 

To explore the effects of the core and crust microphysics 
on the carbon ignition, we choose a core model with M = 
1.6M and z = 0.35, which implies that R = 10.5 km 



and p = 6.2p n 



1.7 x 10 gem 3 . Here p nuc /'«u 



0.1 6 fm~ 3 is the fiducia l nucle ar density. Again follow- 
ing lYakovlev & Haensel (120031) . we use two cases of the 
core neutrino emissivity. The first case is pe v = 3.0 x 
10 19 ergss -1 cm~ 3 (7yi0 9 K) 8 throughout the core. This cor- 
responds to a modified Urea process somewhat suppressed by 
superfluidity. For the second case, we add an enhanced emis- 
sivity pe v = 10 26 ergss- 1 cm- 3 (r/10 9 K) 6 for p > 5.0p nuc . 
This models the operation of an enhanced neutrino emissiv- 
ity from, e.g., a direct Urea process above a density thresh- 
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FIG. 1. — The emergent luminosity from the deep crust, scaled to Lq = 
1 MeV x (M/m n ) for M 17 = 0.1, 0.3, 0.6, and 10.0 (see text). This is for a 
crust with low thermal conductivity. 



old (lLattimer et al These choices change the core 

temperature by an order of magnitude for a given accre- 
tion rate: the proper temperature at the base of the crust is 
r b w 4.6 x 10 8 K(M/M E dd) 1/8 for the first case and T h w 
3.9 x 10 7 K(M/M Edd ) 1/6 for the second. The core temper- 
ature is insensitive to the exact value of the density threshold. 

In addition to the neutrino emissivity of the core, we also 
choose two cases for the crust thermal conductivity. Cal- 
culations of the ionic state of the crust usually presume a 
single-species lattice, with possible impurities. The actual 
crust is formed, however, from the processing of rp-process 
ashes, and it is not known if the ions form an ordered lat- 
tice. We therefore pick two cases for comparison. First, 
that the lattice is disordered; we estimate the electron-ion re- 
laxation time T by s etting the structure factor to unity (see 
lltoh & Kohvamall993l) . In this case T _1 is independent of T 
and proportional to Z 2 /A. For the second case, we compute 
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r e,ph 



r ' = r, + T e j^p, where T e i is the frequency of electron- 



phonon scattering a nd T^jL- is th e frequency of electron- 
impurity scattering (Pote khin et al.lH999l) . for which we set 
Q = (Z 2 — (Z) 2 ) = 100. This is typical of the rp-process 
ashes, but it is larger than that of m atter burned in a super- 
burst (Schatz, Bildsten, & Cumming 2003) or formed at the 
birth of the neutron star (Jones 2004). 

The thermal profiles for these four cases are shown in Fig. [2] 
for M = 3 x 10 17 gs~' (approximately 4 Eddington) and 

Fig- Q] for M = 10 18 gs _1 (approximately Eddington). The 
solutions with the higher core temperatures are for modified 
Urea processes; the ones with lower core temperatures for di- 
rect Urea processes. For each core temperature there are two 
profiles, one for a disordered, amorphous crust {solid lines) 
and one for an impure lattice (dotted lines). 

The squares mark the point at which we estimate un- 
stable ignitio n to occur. To defin e ignition, we use the 
criteria from Cummina & Bildstenj (120011) and set £1212 = 
^pKT/y 2 , where £12,12 is the mass-specific heating rate 
from the reaction 12 C + 12 C, £ CO oi = pKT /y 2 is an es- 
timate of the cooling rate, and | 2/26 is the ratio 
(dj In e coo \) / (dj In 812,12)- This simple approach gives us a 
rough estimate of where ignition occurs and is known for 
X-ray bursts to give results that are roughly in line with nu- 



merical s imulations (IWooslev et alJl2004l) and linear stability 
analyses (Naravan & Hevl 200 3). For the 12 C + 12 C rate, we 
use the formula of Caughlan & Fowler ( 1988) and incorporate 
strong screening according to Ogata, Ichi maru. & van Horrt 
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FIG. 2. — Crust temperatures forM = 3.0 x 10 gs 1 . The squares indicate 
where a mixture with Z( I_ C) =0.1 would unstably ignite. The solid lines are 
for the crust thermal conductivity dominated by scattering; the dotted lines 
indicate a conductivity that is a mix of phonon and impurity scattering. The 
solutions with higher core temperatures correspond include only modified 
Urea processes; the solutions with a cold core temperature include in addition 
an enhanced neutrino emission. 

ForM = 3.0 x 10 17 gs -1 , the recurrence times for 12 C + 12 C 
ignition are, from the hottest profile to the coldest, 1.4, 3.0, 
5.1, and 28 yr with total burst energies of 7.5 x 10 41 , 1.6 x 
10 42 , 2.7 x 10 42 and 1.5 x 10 43 ergss" 1 , respectively. Like- 
wise, forM = 10 18 gs -1 , the recurrence times are 0.08, 0.10, 
0.30, and 1.4 yr and the total burst energies are 1.5 x 10 41 , 
1.7 x 10 41 , 5.3 x 10 41 , and 2.5 x 10 42 ergss" 1 . Neutron stars 
with enhanced neutrino cooling are not compatible with ob- 
served superburst recurrence times and energetics unless the 
crust is amorphous. Note that if the inner crust were to 
form a lattice with a low impurity concentration, Q ~ 1, 
then the 12 C + 12 C ignit ion approaches the case discussed by 
iTaam & Picklum ( 1978), with recurrence times of 10-100 yr. 
For the case of a disordered lattice (solid curves), the ignition 
conditions become nearly independent of the core tempera- 
ture at higher accretion rates. 

4. CONCLUSIONS 

Using parameterized models of the NS core, we have com- 
puted ignition conditions for 12 C that solve for the flux emer- 
gent from the inner crust. There are two principal conclusions. 
First, the flux in the carbon-enriched layer scales roughly as 
(A)/(Z 2 ), so that the superburst ignition is only weakly de- 
pendent on the c omposition of that layer. Sec ond, we confirm 
the scenario of jCumming & Bildsten (2001), but only if the 
crust thermal conductivity is low. In particular, the calculated 
recurrence times are too long if there is both enhanced core 
neutrino emission and the crust thermal conductivity is set 
by phonon and impurity scattering. This is relevant not only 
to studies of neutron star cooling, but also to studies of the 
magn etic field evolution of the cru st (for a recent study, see 
Cumming. Arras. & Zweibel 2004). 

In addition to the recurrence time and burst energetics, 
the cooling of the NS atmosphere, as evidenced by the 
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FIG. 3.— Same as|2] but for M = 10 18 gs~ 



quenching of subsequent type I X-ray bursts, can be used 
to infer y^n, the column de pth of the superburst ignition 
( Cummins & Macbeth 2004). The ignition curve, as esti- 
mated in § |3] gives a unique temperature 7i gn (jign ) ■ Using 
this information thus reduces the uncertainty from the poorly 
constrained mass accretion rate. From Fig. [2] y; gn differs by 
a factor of 3-4- between different crust models for a core with 
only modified Urea processes. 



Perhaps the most intriguing sup erburst source is 
KS 1731-260 (iKuulkers et alJ 12002). for which there 
are independent estimates of the interior temperature and 
crust thermal conductivity. This source was observed to be 
persistent for 12 yr unt il it suddenly w e nt into quiescence 
(Wiin ands et alJ 1200 ll) . iRutledge et aU J2001I) noted that 
because of its long outburst, the crust would be driven 
out of thermal equilibrium with the core, and therefore 
follow-up observations of the quiescent X-ray luminosity 
would measu re the thermal r elaxation of the crust. Such 
observations (Wiinands et al. 2002) revealed a rapid and 
strong decline of the X-ray luminosity. If the quiescent 
emission is solely powered by the cooling crust, this sharp 
decrease in the quiescent luminosity suggests that the crust 
has a high thermal conductivity, consistent with being a 
locally pure lattice, and that the neutrino emission from the 
core is stronger than modified Urea. The calculations in 
this paper, being steady-state, are not directly applicable 
to KS 1731—260, but do indicate that the low quiescent 
emission, combined with the occurrence of a superburst, are 
a challenge to our understanding of accreting neutron stars. 



We thank Andrew Cumming, Bob Rutledge, Hendrik 
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